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ABSTRACT 

I We present the first results from a 1.1mm confusion- limited map of the Great Ob- 

servatories Origins Deep Survey-South (GOODS-S) taken with the AzTEC camera 
on the Atacama Submillimeter Telescope Experiment. We imaged a 270 arcmin^ field 
to a Icr depth of 0.48 — 0.73 mJy/beam, making this one of the deepest blank-field 
surveys at mm-wavelengths ever achieved. Although by traditional standards our 
k>( \ GOODS-S map is extremely confused due to a sea of faint underlying sources, we 

. demonstrate through simulations that our source identification and number counts 

' analyses are robust, and the techniques discussed in this paper are relevant for other 

deeply confused surveys. We find a total of 41 dusty starburst galaxies with signal 
to noise ratios S/N ^ 3.5 within this uniformly covered region, where only two 
are expected to be false detections, and an additional seven robust source candi- 
dates located in the noisier (la « ImJy/beam) outer region of the map. We derive 
the 1.1mm number counts from this field using two different methods: a fluctua- 
tion or "P(d)" analysis and a semi-Bayesian technique, and find that both meth- 
ods give consistent results. Our data are well-fit by a Schechter function model 
with (S", A^3mjy,a) = (1.30;^;25mJy, 160128 mJy"^deg-^ -2.0). Given the depth 
of this survey, we put the first tight constraints on the 1.1mm number counts at 
5*1.101111 = 0.5 mjy, and we find evidence that the faint-end of the number counts at 
'S'sso/iin ^ 2.0 mJy from various SCUBA surveys towards lensing clusters are biased 
high. In contrast to the 870 |i.m survey of this field with the LABOCA camera, we find 
no apparent under-density of sources compared to previous surveys at 1.1 mm; the es- 
timates of the number counts of SMGs at fiux densities > 1 mJy determined here are 
consistent with those measured from the AzTEC/SHADES survey. Additionally, we 
find a significant number of SMGs not identified in the LABOCA catalogue. We find 
that in contrast to observations at A ^ 500 |J.m, MIPS 24 [a.m sources do not resolve 
the total energy density in the cosmic infrared background at 1.1 mm, demonstrating 
that a population of z > 3 dust-obscured galaxies that are unaccounted for at these 
shorter wavelengths potentially contribute to a large fraction (~ 2/3) of the infrared 
background at 1 . 1 mm. 
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1 INTRODUCTION 

Galaxies selected at submillimetre (sub-mm) and millime- 
tre (mm) wavelengths (hereafter SMGs) comprise a pop- 
ulation of dust-obscured starburst or active galactic nu- 
clei (A GN) host galaxie s at high redshift (z > 1; see re- 
view bv lBlain et aJLlbOOd ). With far-infrared (FIR) luminosi- 
ties I/FiR ^ 10^^ L© , these systems appear to be scaled-up 
analogs to the ultra-luminous i nfrared galaxies (ULlRGs ) 
observed in the local Universe jSanders fc Mirabell Il993 ). 
Their FIR to mm spectral energy distributions (SEDs) are 
characterised by thermal dust emission with temperatures 
of Td ^ 35 - 40 K JChapman ct al."2005':'Kova cs et al.ll2006l : 
iPope et all l2006l : [Coppin ct al. 2008), peaking in the FIR 
at rest-frame A ~ 100 [J.m. Due to the steep rise with fre- 
quency of the SEP on the Rayleigh- Jea ns tail {Si, oc 
iDunne et al.ll2000l : lDunne fc Ealedl200ll '). the FIR peak is in- 
creasingly redshifted into the sub-mm/mm observing bands 
with increasing distance, resulting in a strong negative k- 
correction that roughly cancels the effects of cosmological 
dimming with redshift for observations at A > 500 \im. This 
makes SMGs of a given bolometric luminosity equally de- 
tectable between 1 < 2; < 10 at 1.1mm, assuming that a 
sufficient amount of dust can be built up in these systems 
at such early epochs. Extragalactic surveys at sub-mm/mm 
wavelengths have long taken advantage of this unique prop- 
erty and have detected hundreds of starburst galaxies in 
the early Universe, many of which go undetected in even 
the deepest optical surveys due to extreme dust-obscuration 
(e.g. lYounger et akllioOTl bOOgl ). These include deep, large- 
area surveys towards "blank-fields", i .e. regions devoid 



lengths requir e strong luminos it y evo l ution of IR-bright 
galaxies (e.g IScott et akl |2002|. I2OO6I: iGreve et all l2005l : 



Greve et al.l 2004 20081: Laurent et al. 12005: CooDin et al. 
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Bertoldi et all 2007: Scott et al. 20081: Perera et al. 
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Weifi et al.l l2009l: lAustermann et all I2OI0I') and to- 



wards "biased" fields, such as the enviro nments of high- 
redshift radio galaxies (e.g. IStevens et al.i .20 03) and other 
trace rs of high-redshift proto-clusters (e.g. iTamura et al.l 
I2OO9I ). Several small-area surveys towards known clusters 
at moderate redshifts {z < 0.5) have also been carried out 
in order to detect intrinsically less luminous background 
SMGs via amplification through gravitational lensing (e.g. 
Small et al.lll998l. l2()02l: IChaprnan et al.ll2002l : ICowie et all 



20021 : iKnudsen et al.ll2006l . 120081 ) 



The number counts at sub-mm/mm wavelengths pro- 
vide potentially strong constra i nts o n models of galax; 



alaxy 
20051 : 



evolution (e .g. 'Granato et al? '2004'; I Silva et al.l . 

iBaugh et aU [2005: Rowan-Robinson 200^). The hundreds 
of SMGs detected over the past 12 years in large-area sur- 
veys taken with the Submilli meter Common-User Bolome- 
ter Array (SCUBA, 850 ^m; iHoUand et al] Il999l ) on the 
15-m Jai nes Clerk Ma xwell Telescope (JCMT), MAMBO 
fl.2mm: iKrevsa et al.| [l998) on the Institut de Radio As- 
tronomic Millimetrique (IRAM) 30-m telesco pe, Bolocam 
(1.1mm; IClenn et al.l 1 19981 : iHaig et al.l |2004| ) on the 10- 
m Caltec h Submillimeter Ob servatory (CSO), and AzTEC 
(1.1mm; IWilson et all I2OO8I ') on the JCMT and the 10- 
m Atacama Submillime ter Telescope Experiment (ASTE, 
lEzawa et alll2004 [ 20081 ) have put tight constraints on the 
number counts at Sgso/jm ^ 2mJy (Si.imm ImJy). Re- 
gardless of the precise details of various galaxy evolution 
models, the observed number counts at sub-mm/mm wave- 



ICoppin et aUbOOd : lAustermann et al.ll2010l ). The low angu- 
lar resolution (full width at half maximum FWHM ^ 10") of 
these observations makes them confusion-limited at ~ 1 mjy 
(at both 850 ^m and 1.1mm on the 15-m JCMT), limit- 
ing source detections to only the most luminous (Lfir > 
3 X 10^^ Lq) systems with extreme star- formation rates of 
SFR > 5OOM0/yr. Only surveys towards lensing clusters 
have picked out individual sources with lower luminosities, 
and these currently provide only weak constraints on the 
sub-mm number counts at Ss^Oiim ^ 2 mJy. While detect- 
ing a large number of intrinsically fainter systems and es- 
tablishing the link between SMGs and other high-redshift 
galaxy populations must await the improved resolution of 
the 50-m Large Millimeter Telescope (LMT) and the Ata- 
cama Large Millimeter/Submillimeter Array (ALMA), in- 
formation about their number counts at Si.imm < ImJy 
can be gleaned from wide-area, confusion-limited surveys 
with existing facilities, and these can aid in discriminating 
between various galaxy evolution models. 

In this paper we present a 270arcmin^ 1.1mm sur- 
vey of the Great Observatories Origins Deep Survey-South 
(GOODS-S) field taken with the AzTEC camera on the 
ASTE. This is the deepest survey at mm wavelengths ever 
carried out, achieving a root-mean-square (rms) noise level 
of 1(7 = 0.48 - 0.73 mJy. The GOODS-S field represents 
one of the most widely observed regions of sky, with deep 
multi-wavelength data from a number of ground-based and 
spac e-based facilities . This includes X-ray data from Chan- 
dra (|Luo et al.ll2008l ). optical to ne ar-IR photometry from 
the Hubble Space Telescope {HST; ICiavalisco et all 120041 ') . 
near-I R imaging with IS AAC on the Very Large Tele- 
scope iRetzlaff et al.ll2009l '). Spitzer IRAC (Chary et al. in 
prep.) and MIPS (Dickinson et al. in prep.) imaging in 
the mid-IR, sub-mm imaging at 250, 350, and 500 \im with 
the Balloon-bor ne Large Apertur e Submillimeter Tele- 
scope (BLAST; iDevlin et aiTl2009l 'l and at 870 ^.m with 
the LABOCA camera on the Atacama Pathfinder Experi- 
ment CAPEX: iWeifi et al . 2009), and 1.4 G Hz interferometric 
imag i ng with the Very La rge Array (VLA: [Kellermann et al.l 
l2008l : iMiller et al.l [200^ ). Dedicated spectroscopic foUow- 
up of optical sources in this fiel d has also been und erway 
IVanzella et alll2005l . [ioo^ . I2OO8': 'Pop esso et all 12009). Ad- 
ditionally, planned Herschel observations of the GOODS-S 
field will provide ultra-deep FIR/sub-mm imaging at 100, 
250, 350, and 450 [im. This suite of multi-wavelength data 
is essential for the identification of counterparts to the SMGs 
and for the characterisation of the properties of these galax- 
ies. 

This paper is organised as follows: in Section [5] we de- 
scribe the observations of the GOODS-S field carried out 
with AzTEC on ASTE. In Section [3] we summarise the 
data reduction methods. We present the 1.1mm map and 
source catalogue in Section [l] and we describe simulations 
carried out to characterise the number of false detections, 
survey completeness, and degree of source blending in the 
map in Section [5] We derive the 1.1mm number counts 
from this survey in Section [6] and compare them with the 
number counts determined from SCUBA lensing cluster sur- 
veys and existing blank-field surveys at 1.1 mm wavelengths. 
We present a comparison between the AzTEC/GOODS-S 
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data and the 870 [im data from LABOCA in Section [7] We 
discuss the contribution to the cosmic infrared background 
(CIRB) at 1.1mm from the radio and mid-IR galaxy pop- 
ulations in Section [S] and we close with a summary of our 
results in Section (9] Several upcoming papers involving this 
data-set are underway, including the identification of ra- 
dio and mid-IR counterparts to the SMGs (Scott et al. in 
prep.), a study of the 1.1mm properties of BzK-selected 
galaxies (Welch et al. in prep.), a comparison between the 
AzTEC/ASTE GOODS-S data and the sub- mm maps from 
BLAST (Aguilar et al. in prep.), and a study of the X-ray 
properties of SMGs (Johnson et al. in prep.). 



2 OBSERVATIONS 

We imaged a 26 x 20 arcmin^ field centred at right ascen- 
sion and declination (RA, Dec) = (03''32™30", -27°48'20") 
at 1.1mm using AzTEC on the ASTE. The central 19 x 
14arcmin'^ region, where the coverage is uniform, encom- 
passes the entire GOODS-S region mapped by the Spitzer 
Space Telescope. The observations were carr ied out using the 
N-CO SMOS3 network observation system (jKamazaki et al.l 
bOOSl ) from July 15 to August 6 during the 2007 Chilean win- 
ter under excellent observing conditions, with T220 ~ 0.05 on 
average, and T220 < 0.06 70% of the time (zenith opacity at 
220 GHz reported by the ASTE tau monitor). A total of 
52 hours of observing time excluding pointing and calibra- 
tion overheads was devoted to this survey. During the 2007 
season, 107 (out of 144) of the AzTEC bolometers were sta- 
ble with high sensitivity. The point spread function (PSF) of 
each detector was measured via beam-maps on Uran us, Nep- 
tune, and 3C279 as described in lWilson et al.l l|2008l ) and has 
a FWHM of 30" ± 1" and 31" ±2" in azimuth and elevation, 
respectively (the theoretical beam size is 27" FWHM). The 
full array subtends a circular field-of-view with a diameter 
of 8'. 



2.1 Scan Strategy 

We used a continuous scanning strategy which traces a mod- 
ified Lissajous pattern on the sky: 

ARA = 5.5' • sin(a ■ t + 0.25) + 2.0' ■ sin(a • t/30) 
ADec = 7.5' •sin(6-t) + 2.0' •sin(6-t/30), (1) 

where a/b = 8/9 and ARA and ADec are physical coordi- 
nates relative to the field centre. The actual values of a and 
b were scaled to limit the peak scanning velocity to 300"/s, 
and a rotational angle 20° West of North was used in order 
to align our map with that of the Spitzer IRAC/MIPS cov- 
erage of GOODS-S. A single observation took 42 minutes to 
complete, and we obtained a total of 74 observations of the 
GOODS-S field. 

The benefit of Lissajous scanning, in addition to at- 
taining excellent cross-linking and uniform coverage in the 
map, is that we avoid large telescope accelerations that can 
induce systematics in the detector signals as well as com- 
promise the pointing accuracy of the telescope. Such effects 
are often seen in images taken in raster-scan mode, where 
1/3 — 1/2 of the data taken during times when the telescope 
reverses direction must be discarded (e.g. IScott et al.ll200j : 



IPerera et al.ll2003 'l. Lissajous scanning, on the other hand, 
results in nearly 100% observing efficiency. 

2.2 Pointing Corrections 

We make small corrections to the telescope pointing model 
by routinely observing the bright point source J0455— 462 
(Si.imm ~ 1.5 Jy, variable) every two hours before and after 
each block of GOODS-S observations. We measure pointing 
offsets by fitting the 4 x 4arcmin^ maps of J0455— 462 to 2- 
dimensional Gaussians, and we linearly interpolate these off- 
set corrections temporally and apply them to the GOODS-S 
data. The random pointing error in the final GOODS-S map 
is < 1" (see Section [421). 

2.3 Flux Calibration 

The flux conversion factor (FCF) used to convert the raw de- 
tector signals to flux density units was determined by beam 
map observations on Uranus, Neptu ne, and 3C279 , take n 
1 — 2 times per night as described in IWilson" et al.l l|2008h . 
The flux densities of Uranus and Neptune at 1.1mm were 
calculated from their frequency-depend e nt br ightness tem- 
peratures reported in ICriffin fc OrtonI (|l993h and ranged 
from 43 — 52 Jy and 18 — 20 Jy, respectively, during the 
AzTEC/ASTE 2007 observing season. The flux density of 
3C279, which is highly variable, is adopted from the Sub- 
millimeter Array flux density archivqj and ranged from 
7.0 — 9.4 Jy at 1.1 mm during this time. We remove the re- 
sponsivity factor from the detector signals and correct for 
extinction by modelling both as a linear fun ction of the de- 
modulated DC- level Csee lWilson et al.l l2008''). 

The measured FCF varied significantly from night to 
night, resulting in a la scatter of 17% over the entire ob- 
serving run. We have identified the source of this scatter as 
the changing focal point of the telescope with environment 
temperature: the FCF decreases as the measured FWHM of 
the beam increases. Since real time corrections to the sub- 
reflector position were not possible, we use the same-night 
measurement of the FCF to calibrate each observation. To 
estimate the total calibration uncertainty, we determine the 
standard deviation in the measured flux densities from the 
68 pointing observations of J0455— 462, which is 8%. Since 
this source is known to be variable, this gives a conserva- 
tive upper limit to our calibration uncertainty. Combining 
this in quadrature with the 5% absolute uncertainty on the 
flux densities of Uranus and Neptune (|Griffin fc Ortonlll99^ ) 
gives a total upper limit to the calibration error of 9%. 



3 DATA REDUCTION 

We reduce the 1.1 mm data in a ma nner that is nearly iden- 
tical to that described in detail in IScott et al.l (|2008h . We 
summarise the steps here and note the differences. The raw 
time-stream data, which consist of all bolometer signals and 
pointing data stored as a function of time, is first scanned 
for "spikes" (defined as any > 7a jump between sequential 
detector samples), which are caused by instrumental glitches 

^ http://smal.sma.hawaii.edu/caUist/caUist.html 
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or cosmic ray strikes. These data and nearby samples, which 
amount to < 0.1% of the total time-stream data, are flagged 
and discarded from the data-set. We group the remaining 
samples into 10-second intervals, and then "clean" each 10- 
second group using a principal component analysis (PCA) to 
identify and remove the common-mod e atmospheric signal 
l|Laurent et al.ll2005l : IScott et al.ll2008l ). 

For AzTEC maps taken in ra ster-scan mode (e.g. 
IScott et al] I2OO8I : IPerera et all I2OO8I ) data samples from 
the same individual scan (a single pass of the telescope 
across the sky) were grouped together for PCA clean- 
ing. Since there is no such natural division for continuous 
Liss ajous-scan maps, w e chose a 10-second interval group- 
ing. iPerera et al.l (|2008l ) showed via a statistical correlation 
analysis that this PCA cleaning technique using time inter- 
vals ranging from 5 — 15 seconds provides a good balance 
between using a sufficient number of samples to determine 
the bolometer-bolometer correlations and being on a short 
enough time scale so that the slower electronics related low- 
frequency drifts can be effectively removed. We have verified 
that cleaning on 5 — 20-second intervals gives consistent re- 
sults for AzTEC/GOODS-S. 

After cleaning the time-stream data, the bolometer sig- 
nals are calibrated and binned into 3" x 3" pixels to make a 
map for each separate observation. These 74 maps are then 
co-added and wiener filtered to suppress correlated large- 
scale structure from residual atmosphere and pixel-to-pixel 
variations on scales smaller than the bea m in order to opt i- 
mise the map for point source detection (jScott et al.ll200^ . 
In addition to this filtered map, we track the effects of PCA 
cleaning and filtering on a model point source (referred to 
hereafter as the point source kernel) . The smoothing slightly 
broadens the FWHM of the beam: fitting a 2-dimensional 
Gaussian to the point source kernel results in a FWIfM of 
34.6" and 34.3" in RA and Dec, respectively. We also gener- 
ate 100 noise maps, each a realisation of the signal-free noise 
in the GOODS-S map, by "jackknifi ng" the time - stream sig- 
nals in the same manner described in lScott et all l|2008l ). The 
point source kernel and noise realisations are used later in 
the analysis for simulation purposes (Section[S]l and number 
counts determination (Section |6}. 



4 1.1 MM MAP AND SOURCE CATALOGUE 
4.1 AzTEC 1.1mm Map 

The fuU 512arcmin^ 1.1mm map of GOODS-S taken with 
AzTEC on ASTE is shown in Figure [1] Since we must con- 
sider only a region with uniform noise for most of the anal- 
ysis presented in this paper, we define a 270 arcmin'^ "uni- 
form coverage region" , where the coverage is greater than or 
equal to 50% of the maximum coverage in the map. The la 
rms noise in this region ranges from 0.48 — 0.73 mjy/beam, 
making this the deepest contiguous region ever mapped at 
1.1 mm. The scanning strategy that we used to map this field 
results in two separate patches where the coverage is deepest 
(0.48 mjy /beam), whereas the center of the map is slightly 
shallower (0.56 mJy /beam, see rms contours on Figure O. 
The noise determined from the jackknifed noise realizations 
is Gaussian distributed with la = 0.57 mJy/beam. 



4.2 Astrometry 

To check for a residual astrometric offset after applying 
the pointing model (see Section |2.2|I we stack the 1.1mm 
GOODS-S map at the positions of radio sources in this field, 
whose positions are known to much better than 1". From the 
VLA 1.4GH[z sur vey of the Extend ed Chandra Deep Field- 
South (ECDF-S; iMiller et all 120081 ). which reaches an rms 
noise level of 7 — 10 /iJy, we extract a ^ 4cr radio source 
catalogufl and stack the AzTEC map at the positions of 
the 219 radio sources that lie within the 50% uniform cov- 
erage region. The resulting stacked AzTEC map shows a 
clear peak with signal to noise S/N — 8, and its centroid 
is offset by -6.3" ± 2.1" in RA and -3.9" ± 2.1" in Dec. 
We have verified this result by stacking at the locations 
of 1185 MIPS 24 \xm sources detected with S/N > 10 in 
the Spitzer GOODS-S survey (Dickinson et al. in prep.), 
which gives a S/N = 11 detection with a centroid offset of 
-6.3" ± 1.6" in RA and -0.7" ± 1.6" in Dec, consistent with 
the radio stacking results. We thus apply an astrometric cor- 
rection of (ARA,ADec) = (6.3", 0.7") to the AzTEC map 
and 1.1mm source positions (we favor the offsets measured 
from the 24 \im stacking due to the higher S/N of the peak 
detection) . 

This systematic offset represents the average pointing 
offset between J0455-462 and the GOODS-S field over 74 
observations. The scatter in this offset, or the pointing jitter, 
will manifest itself as a broadening of the source in excess 
of the point source kernel. We estimate the random point- 
ing error in the AzTEC map from the brightest AzTEC 
source in this field, AzTEC/GSl. We adopt a simple model 
which consists of the convolution of the ideal point source 
kernel with a 2-dimensional Gaussian with standard devi- 
ation (cr RA, cTDec), where ctra and ffp cc are the la ran dom 
pointing errors in RA and Dec (see IScott et al.l |2008| . for 
a full description of this measurement). The high S/N of 
AzTEC/GSl (11. 60", see Table [l| allows a clean measure- 
ment and provides a strong constraint on the random point- 
ing error, since the signal from this single source can only 
be broadened with respect to the point source kernel due to 
the scatter in the pointing model. The maximum likelihood 
estimate for the random pointing error from AzTEC/GSl 
is (o"RA,o"Doc) ~ (0.5", 0.1"); however, the distribution of 
(c7RA,o"Dec) is Very flat out to 1", then falls off steadily. 
From this we conclude that the random pointing error in 
the AzTEC/GOODS-S map is < 1". 

While the stacked 1.1 mm signal for the radio and MIPS 
24 (i.m populations is valuable for determining the astromet- 
ric offset, we find that the stacked signal cannot be used for 
determining the pointing jitter as it is significantly broader 
than the point source kernel, with la « 6.8" and 15" for the 
stack on the radio and 24 ^m sources, respectively. This ad- 
ditional broadening is likely caused by other effects such as 
confused/blended sources in the 1.1 mm image or clustering 
of the radio and 24 ^.m sources. 



^ The radio source catalogue is produced using the SAD pro- 
gram in the Astronomical Image Processing System (AIPS; 
http://www.aips.nrao.edu/), which uses 2-dimensional Gaussian 
modelling to identify sources and produce photometry results. 
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Figure 1. The AzTEC 1.1mm map of the GOODS-S field. The sohd contour shows the boundary of the 270arcmin'^ 50% uniform 
coverage region to which most of the analysis of this field is restricted. The dashed contours (innermost to outermost) indicate noise 
rms levels of 0.51 and 0.57 mjy/beam. The circles located within the solid contour (diameter = 2 X FWHM of AzTEC on ASTE = 60") 
indicate the positions of ^ 3.5cr source candidates, labeled in order of decreasing S/N of the detections. The circles located outside of 
the solid contour indicate robust > A.5cr sources within the noiser regions of the map. See the on-line journal for a colour version of this 
Figure. 



4.3 Source Catalogue 

We identify point sources in the 1.1 mm S/N map by search- 
ing for local maxima within 15" of pixels with S/N J5 3.5 
(see lScott et~al]|2008l , for a more detailed description). The 
40 source candidates located within the uniform coverage 
region that meet this criterion are listed in Table [1] in or- 
der of decreasing S/N of the detection. We note that the 
number of selected source candidates is independent of the 
window size for grouping high S/N pixels, for values rang- 
ing from 3"to 45". Table [T] includes both the 1.1mm flux 
densities and la errors measured from the map, as well 
as the bias-corrected flux densiti es estimated using a semi- 
Bayesian technique (Section fiOl ICoppin et all |2005| . l2006l : 
lAustermann et aDl2009l . |2010( ). This flux-bias correction ac- 
counts for the fact that the measured flux densities of mm- 
selected galaxies, which are generally detected at low S/N, 



are preferentially "boosted" due to the steep luminosit y dis- 
tribution of the population (e.g. iHogg fc Turneij|l998l ). 

We flnd an additional seven robust source candidates 
that are located outside of the uniform coverage region, but 
are detected with high S/N (^ 4.5). These sources are listed 
in Table [T] for the beneflt of future studies; however, we do 
not address the general properties of these sources located 
outside of the uniform coverage region (e.g. number of false 
detections, completeness) nor include them in our number 
counts estimation in Section |6] 

Several of the source candidates appear extended in 
the 1.1mm map, most notably, AzTEC/GS2 (labelled as 
2.1 and 2.2 in Figure [!}. To separate the components of 
AzTEC/GS2, we fit the 1.1mm map in the neighbourhood 
of this source to a 2-component model, where each compo- 
nent is a scaled version of the point source kernel. The best- 
fit positions and flux densities are listed in Table [1] Given 
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the comparatively modest S/N of the other AzTEC sources 
which are potentially extended, we cannot fit a model to 
these sources in order to separate them into multiple com- 
ponents. 



5 CHARACTERISATION OF A DEEPLY 
CONFUSED MAP 

Source confusion arising from a significant underlying pop- 
ulation of faint sources that are individually undetected 
in a given sur vey is a fu nction of source density and 
resolution (e.g. iHoggj l200ll . and references therein), and 
is important to consider for our AzTEC/GOODS-S sur- 
vey. Using the standard rule of thumb (one source per 
30 beams jj and the blank-fiel d 1.1mm source counts from 
the AzTEC/SHADES survey l|Austermann et al.llioiol '). the 
confusion limit given the 30" FWHM ASTE beam is 
2.0 mjy. Since our AzTEC/GOODS-S survey is 3 - 4 times 
deeper than the formal confusion limit, we use specialised 
simulations to characterise the effects of source confusion on 
the properties of our GOODS-S map, including the number 
of false detections, the survey completeness, the positional 
uncertainty distribution, and the degree of source blending. 

For these purposes, we generate sky realisations of the 
GOODS-S field using the noise maps and simulated point 
sources, each modelled as the point source kernel scaled by 
the source flux density. These simulated galaxies are in- 
jected at random positions drawn from a uniform distri- 
bution into the noise maps, where their flux density dis- 
tribution is described by the best-fit Schechter function to 
the AzTEC/GOODS-S number counts (Section [6)1 trun- 
cated at 5i.imm = 0.1 mJy. These maps, referred to here- 
after as "fully simulated maps", provide a realistic model 
of the mm-galaxy population in the GOODS-S field as ob- 
served by AzTEC by properly including the effects from our 
data reduction methods. We use these fully simulated maps 
throughout this paper to investigate various properties of 
our map. 

One might ask if the assumption of a spatially uni- 
form population of sources in our simulations is justified. 
To date there are only weak constraints on the clustering 
properties of bright SMGs derived from incomplete and/or 
small-number samples , with significant disagreement among 
various surveys (e .g. Isiain et all l2004l : iGreve et all l2004l : 
IScott et al.l I2OO6I : IWeifi et al.fbOOgl '). Furthermore, bright 
sources like those probed in these surveys have low number 
density and so will not bias the map characteristics probed 
in our simulations. Rather, it is the faint, high-density, un- 
explored population of SMGs that can bias our estimates. 
Unfortunately, until better resolution for sub-mm/mm sur- 
veys is available (for example AzTEC on the LMT), we 
have no sufficient observational evidence to guide us. While 
these faint SMGs may have similar clustering properties as 
other high-redshift galaxy populations (e.g. LBGs, BzKs, 



^ B y deriving a general form alisation of source confusion statis- 
tics, l^^uchT&JshO 1I2OO4I ') demonstrate the basis for the "one 
sourc e per 30 beams" rule of thumb for estimating the confusion 
limit. iTakeuchi fc Ishiil \2m4 ) also show that the actual confusion 
limit is strongly dependant on the steepness of the source counts, 
rather than simply the source density. 



EROs) it is likely that the projected 2-dimensional cluster- 
ing strength of SMGs is very different from optically selected 
high-redshift galaxies, since sub-mm/mm surveys sample a 
much larger volume of space. In the absence of observational 
constraints on the clustering of faint SMGs we choose to 
model the source population with a uniform spatial distri- 
bution. 



5.1 False Detections 

Given the modest S/N of the source candidates, we expect 
some fraction of the AzTEC sources in GOODS-S to be spu- 
rious. One method commonly used to estimate the number 
of false positives is to run the source-finding algorithm on 
the negative of the signal map. However, the PGA cleaning 
used for atmosphere removal leaves the signal map (and the 
point source kernel) with a mean of zero. Every real source 
that increases the number of positive pixels in the map will 
thus also increase the number of negative pixels, so the num- 
ber of false detections measured from the negative of the real 
map will be overestimated. Instead, we estimate the number 
of false detections by identifying the number of "sources" ex- 
tracted from the 100 pure noise realisations. The number of 
false detections expected as a function of limiting S/N ratio 
is shown in Figure [5] (solid curve, diamonds). At ^ 3.5(t, we 
expect ~ 2/40 sources in our catalogue (5%) to be spurious. 
Above ^ 4.25(T, the number of AzTEC sources expected to 
be false positives is consistent with zero. 

This estimate however provides only an upper limit to 
the number of spurious detections. In the real map, the neg- 
ative bias in the pixel fiux distribution from the addition 
of sources decreases the number of high-significance positive 
noise peaks in the map. This effect was first demonst rated 
for the AzTEC/GOODS-N survey (|Perera et al.ll2008 ll and 
is even more pronounced for our confusion-limited GOODS- 
S map. To demonstrate this effect, we run the source-finding 
algorithm on 600 fully simulated maps: for each detected 
source, we search within a 17" radius to identify the cor- 
responding input source, and detected sources that cannot 
be traced back to an input source are false positives. We 
choose a 17" search radius to ensure recovery of an input 
source, since a source detected at S/N ^ 3.5 has a proba- 
bility < 0.1% of be ing detected >17" from its true position 
jlvison et al.l|20 (y?). However given the high source density 
at faint fiux densities, there is a non-negligible probability 
that an intrinsically faint source will be located within 17" 
of a random position in the map, and so this calculation will 
underestimate the number of false positives. We account for 
the fraction of random associations as follows: 1) for each 
of the 600 fully simulated maps, we select 2, 000 random 
positions within the 50% uniform coverage region, identify 
input sources located within 17" of these random positions, 
and from this determine the probability P(> 17"|5'i) that a 
source with intrinsic fiux density 5"; will be located within 
17" of a random position in the map; 2) for each detected 
source in the simulated maps that can be traced back to 
an input source with intrinsic fiux density 5*;, we generate 
a random number, Ptcst, between and 1 from a uniform 
distribution, and we consider the output-input source pair 
a false association (and hence the output source a false pos- 
itive) ff Ptcst < P(> 17"|S0. 

The number of false positives estimated from these fully 
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Table 1. AzTEC/GOODS-S source candidates. The columns give: 1) AzTEC source identification; 2) source name; 3) S/N of the 
detection in the AzTEC map; 4) measured 1.1mm flux density and error; and 5) deboosted 1.1mm flux density and 68% confidence 
interval (Section |62J. The sources above the horizontal line are detected with S/N 3.5 writhin the uniform coverage region of the 
AzTEC/GOODS-S map, where two are expected to be false positives. The S/N 4.5 sources listed below the horizontal line are located 
outside the uniform coverage region. 
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simulated maps as a function of limiting S/N is shown in 
Figure [2] (dashed curve, squares). From this estimate we ex- 
pect at most one of the 40 ^ 3.5(t sources in our catalogue 
to be spurious. At ^ 3. Oct, the number of false positives esti- 
mated from these fully simulated maps is significantly lower 
than that estimated from pure noise realisations (1 — 2 ver- 
sus ^ 7), suggesting that we can comfortably extend our 
source catalogue to lower S/N detections. Ifowever, these 
simulations do not include the effects of source clustering. 
If the mm-galaxy population is strongly clustered on small 
angular scales, the strength of the negative bias in the pixel 
flux distribution would vary from region to region due to the 
variance in the source density, and thus the number of posi- 
tive noise peaks (i.e. false positives) would be lower (higher) 
in the more (less) densely populated regions of the map. 
Since the clustering properties of the mm-galaxy population 
(including intrinsically faint sources) are not well known, we 
prefer to quote the values determined from the pure noise 
realisations as a conservative estimate of the number of false 
detections expected in our catalogue. 

5.2 Completeness 

The detection rate for a given source flux density is af- 
fected by both Gaussian random noise in the map and con- 
fusion noise from the underlying bed of faint sources. To 
account for both effects, we estimate the survey complete- 
ness by measuring the recovery rate of point sources with 
known flu x densities inserte d into the real signal map, as de- 
scribed in lScott et al.l (|2008l ). For flux densities ranging from 
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Figure 2. The expected number of false detections in the 
AzTEC/GOODS-S map as a function of limiting S/N. The solid 
curve and diamonds show the number of false detections esti- 
mated from the number of peaks detected in pure noise reali- 
sations with significance ^ S/N. The dashed curve and squares 
indicate the expected number of false positives determined from 
fully simulated maps. The pure noise realisations provide only 
a conservative upper limit to the number of false detections ex- 
pected. 



0.1 — 8.0 mjy, we input 2000 sources per flux density one at 
a time into the GOODS-S map, each time randomly select- 
ing the source position. Using the standard source-finding 
algorithm, an input source is considered recovered if it is 
detected in the map within 17" of its input position with a 
significance of ^ 3.5a. We exclude samples where the simu- 
lated source was input or extracted less than 17" from a real 
J5 3.5(7 source. The survey completeness is shown in Figure 
[3] (data-points with binomial error bars). The survey is 50% 
complete at 2.1 mJy, and 95% complete at 3.5 mJy. 

The sea of faint sources below the detection threshold 
adds confusion noise to the AzTEC/GOODS-S map. This 
additional noise reduces the map's sensitivity to individual 
sources and its survey completeness over a range of fiux den- 
sities. An accurate estimate of the completeness is essential 
for correcting the observed number counts for this field. In 
the standard Bayesi an method for extracting numbe r counts 
from AzTEC maps (jAustermann et aLllioog . I2OI0I ). survey 
completeness is estimated by injecting sources of known flux 
density one at a time into noise realisations and determining 
their recovery rate. While this method does not take confu- 
sion noise into account, this estimate was found to be con- 
sistent with that measured from the real signal map using 
the method described in the previous paragraph for several 
AzTEC surveys on the JCMT, where the angular resolution 
is better (FWHM = 18" on the JCMT versus 30" on the 
ASTE) and the la depths of the surveys are 1.0 mJy (e.g 
IPerera et all I2OO8I : lAustermann et all |2009| . I2OI0I ) , demon- 
strating that confusion noise was not signiflcant for these 




Flux Density (mJy) 

Figure 3. The survey completeness for S/N ^ 3.5 AzTEC 
sources in GOODS-S. The data-points and 68% confidence bi- 
nomial error bars show the completeness estimated by inserting 
sources of known flux density one at a time into the real signal 
map. The solid curve shows the completeness estimated by fully 
simulated maps. The dashed curve shows the completeness esti- 
mated by inserting sources of known flux density one at a time 
into pure noise realizations. This latter estimate does not account 
for the effects of confusion noise and hence results in an overes- 
timate of the survey completeness for the range of flux densities 
shown here. 
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observations. In contrast, we find tiiat tlie completeness 
estimated from noise-only maps significantly over-predicts 
the survey completeness for our deeper, confusion-limited 
GOODS-S map (Figure [S] dashed curve). 

To verify that this difference arises from confusion 
noise, we next estimate the survey completeness from 10, 000 
fully simulated maps. For each ^ 3.5(t source detected in 
these simulated maps, we identify the brightest input source 
within 17" of the output source position. We bin all detected 
sources by their input flux densities, and the completeness is 
calculated by the ratio of the number of recovered sources to 
the total number of input sources per flux bin. The complete- 
ness estimated from these fully simulated maps is shown 
as the solid curve in Figure [S] This estimate agrees quite 
well with that from the single-input source simulations us- 
ing the real GOODS-S map. The discrepancy between the 
two methods at Si.imm ^ 1.5 mjy likely arises from small 
imperfections in the assumptions we use to identify input- 
output pairs. For example, the single-input source simula- 
tions may slightly overestimate the completeness at low flux 
densities due to cases when the input source is inserted close 
to (but > 17") a bright mm-galaxy in the real map. Still, 
despite the very different methods used for these two differ- 
ent completeness estimates, they differ by ^ 4% at all flux 
densities. 



5.3 Positional Uncertainty 

The large beam combined with the low S/N of the detections 
results in a large positional error on the locations of sub- 
mm/mm-detected sources due to the effects of random and 
confusion noise in the map. We use the simulations described 
in l5.2l where a single source of known flux density is inserted 
into the GOODS-S map one at a time, to determine the 
distribution of input to output source distances as a function 
of detected S/N. The probability P(> 0- S/N) that a source 
will be detected outside of a radial distance 6 of its true 
position is shown in Figure |4] for three sample S/N bins. 
For compari s on, t he analytical solutions determined from 
llvison et al.l (|2007l ') are shown as the solid (3.5 < S/N < 
3.75), dashed (4.5 < S/N < 4.75), and dotted (5.5 < S/N < 
5.75) curves, assuming that the FWHM of the AzTEC beam 
is 34" (i.e. the width of best-fit Gaussian to the filtered point 
source kernel). The analytical and empirical distributions 
show broadly the same trend. 

5.4 Source Blending 

Given the depth and low angular resolution of the 
AzTEC/GOODS-S survey, some fraction of the ^ 3.5cr 
sources in the map are expected to be the combined signal 
from two or more individual galaxies blended together. To 
estimate the fraction of the sources in Table[T]that are actu- 
ally the blend of ^ 2 individual galaxies, we take the ^ 3.5a 
sources detected from 600 fully simulated maps and trace 
each one back to all input sources located within 17" of the 
output source position. The fraction of detected sources that 
cannot be traced back to any input source is 0.8%: these rep- 
resent the "false positives" estimated from fully simulated 
maps as described in Section [5. II Sources that map back to 
only one input source are considered "single sources" , while 
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Figure 4. The positional uncertainty distribution for 
AzTEC/GOODS-S source candidates. The data-points and 
error bars show the probability P(> B; S/N) that an AzTEC 
source detected with a given S/N ratio will be found outside a 
radial distance 6 from its true location as determined from sim- 
ulation (Section 15.311 . T he cu rves show the analytical expression 
derived in llvison et all ll2007ll : solid, 3.5 < S/N < 3.75; dashed, 
4.5 < S/N < 4.75; and dotted, 5.5 < S/N < 5.75. The empirical 
and analytical distributions show broadly the same trend. 



those that can be traced back to two or more input sources 
are considered "blended sources" . 

With this simple definition, we would expect nearly all 
sources to be blended given the large beam and the high 
source density of SMGs. However, a very faint source nearby 
a relatively bright source (for example, a 0.2 mJy source lo- 
cated 10" from a 3.0 mJy source) contributes a negligible 
amount to the summed signal. We want to avoid counting 
cases like these - where the brighter of the two sources com- 
pletely dominates the total signal - as a blended source. 
As a more practical definition, we only consider a pair of 
nearby sources to be a blend if the contribution from each 
source to the summed signal is comparable. Using the input 
source fiux densities and relative separations for all simu- 
lated sources within 17" of an output source, we model the 
beam-smoothed noiseless signal from the sum of these point 
sources and measure the peak fiux density. If the fractional 
contribution to the summed fiux density for an individual in- 
put source at the location of the peak is ^ 70%, we consider 
this a single source; otherwise, the detection is considered 
a blended source. With this definition, 25% (10/40) of the 
AzTEC sources listed in Table [T] are expected to be ^ 2 in- 
dividual galaxies blended by the large beam. This fraction 
reduces to 15% if we define a detection as a single source 
when it contributes ^ 60% to the summed noiseless fiux 
density. These results represent lower limits to the fraction 
of blended sources, since this fraction would be even higher if 
the SMC population is significantly clustered. We note that 
by design, this fraction does not depend on the limiting S/N 
threshold (at least for S/N > 3.5). 
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5.5 Summary of Confusion Effects 

We summarise our primary findings on the effects of confu- 
sion on the map properties: 

• When the signal from confused sources becomes com- 
parable to that from random noise fluctuations, the number 
of pure-noise positive peaks in the map decreases. Thus the 
number of false positives as a function of limiting S/N ra- 
tio is considerably lower than that estimated from peaks in 
noise-only realisations. 

• In terms of the survey completeness for a given cata- 
logue, the confused signal acts as an additional noise compo- 
nent that reduces the map's sensitivity to individual galaxies 
over a wide range in intrinsic source flux density. It is impor- 
tant to include the effects of confusion when estimating the 
survey completeness for use in number counts calculations 
(Section Ell). 

• A potentially large fraction of sources detected in a 
confusion-limited map are actually multiple galaxies blended 
by the large beam. Care must be taken to understand the 
possible biases this could cause in number counts estima- 
tions (Section 16. 2p . 



6 NUMBER COUNTS 
6.1 Fluctuation Analysis 

Due to the exceptional depth reached by this survey, the 
mm-emission from faint SMGs has a striking effect on the 
flux density distribution in the map. As discussed already, 
the method used to remove low-frequency modes leaves the 
mean of the map and the point source kernel equal to zero, 
and every mm-source adds both positive and negative flux 
to the map. This is demonstrated in the left panel of Fig- 
ure O which shows the histogram of flux density values in 
the AzTEC/GOODS-S map (with Poisson error bars). The 
dashed curve shows the distribution of pixel fluxes averaged 
over 100 jackknifed noise realisations of this fleld, and is very 
Gaussian. The flux distribution in the real map on the other 
hand is skewed by the presence of mm-sources. While this 
effect makes the identiflcation of individual galaxies chal- 
lenging, we can use the distribution of flux values in the 
map to perform a fluctuation analysis, or more commonly 
referred to as a "P(d)" analysis, in order to determine the 
number counts distribution for this field. As this technique 
is independent of the level of confusion in the map (at the 
expense of being somewhat model dependent), we can use 
it on the AzTEC/GOODS-S map to provide strong con- 
straints on the SMG number counts at faint fiux densities 
(5*1. imm < ImJy) well below the 1.1mm confusion-limit. 

The fluctuation analysis is carried out as follows: using 
a parametrised model of the number counts, we generate 100 
simulated maps as described in Section [5] and compare the 
flux density distribution averaged over these simulated maps 
to that of the real GOODS-S map. For this single model, we 
calculate the comparison metric: 



ln(£) = ^ rrii ~ di + di ■ ln{di/mi) 



(2) 



corresponding quantity for the GOODS-S map. This pro- 
cess is repeated over a grid in parameter space to find the 
minimum of the above metric, which occurs at the best-fit 
model. Minimising this metric is equivalent to finding the 
maximum likelihood for the case that all histogram bins 
follow independent Poisson distributions. Note that we do 
not attempt to model effects of source clustering on the 
pixel flux distributions, since the clustering properties of the 
SMG population are not well known. This method is simi- 
lar iri_£rirra£l«to_tl^ frequentist approach used 
bv lPerera et al.l l|2008h to determine number counts for the 
AzTEC/GOODS-N survey; however, here we consider the 
full flux density histogram in order to extract information 
about the faint source population. 

For this analysis we choose a Schechter function model 
given by: 



ds = ^^-^^ [s^ 



a + l 



exp 



-(g-SmJy) 
S" 



(3) 



where 



djv 
"3? 



is the differential number counts as a function 



of intrinsic 1.1mm flux density S, and {S',Nsnijy,a) are 
the free parameters. While there are many forms of the 
Schechter function published in the literature, we prefer to 
fit to A'smjy because it reduces the degeneracies between the 
parameters, and it has been used in the n umber c ounts anal- 
yses of several previous AzTE C surveys (|Perera et al.. 200^ : 
lAustermann et"al ] |2009l . l2010[ l. making it straightforward to 
compare the results. We prefer a Schechter function model 
over that of a single power-law because it allows for a nat- 
ural steepening of the counts at high fiux densities, which 
has been co nfirmed by large-area surveys at sub-mm/m m 
wavelengths (|Coppin et al. I I2OO6I : lAustermann et ai]|2010l ') - 
though we are unlikely to see this steepening in the GOODS- 
S number counts given the small survey area. Since our data 
cannot provide strong constraints on a 3-parameter fit, we 
fix q: = —2: a value that is consistent with estimates from 
previous surveys f e.g.lCoppin et al.ll2006l : [ Perera et al ] |2008l : 



where rrii represents the average number of pixels in the 
i*'^ flux density bin from the model and di represents the 



lAustermann et al.ll2009l . I2OI0I ) 

Since the Schechter function increases to infinity as S 
goes to zero, we must assume some minimum fiux den- 
sity cutoff, Smin, for the population. A practical minimum 
flux limit is imposed by the data itself: at the flux den- 
sity corresponding to where the number density of sources 
is ~ 1 per beam, adding fainter sources will not alter the 
flux density distribution in the map. Assuming the best- 
fit m odel to the AzTEC/SHADES data l|Austermann et all 
l2010l . which currently provides the tightest constraints on 
the blank-field 1.1mm number counts), Smin ~ 0.1 mjy. 
While it is not known whether the number density of sources 
for the SMG population turns over and starts to decrease 
somewhere below 1 mJy, the counts at the faint end of the 
850 \im SCUBA galaxy population dete r mined from lens- 
ing cluster surveys (e.g. j cpwie et al.ll2002l : ISmail et al.ll2002l : 
iKnudsen eralll2006l . I 20od ) continue to rise out to Ss^Ofim ^ 
0.2 mJy, giving some reassurance that a 1.1 mm fiux cutoff of 
Smin ~ 0.1 mJy is reasonable. We use ^min = 0.1 mJy in gen- 
erating all simulated maps discussed in this paper; however, 
we have tested values ranging from Smin = 0.05 — 0.3 mJy 
and have found that this does not affect the results from the 
fluctuation analysis. 

Using 0.1 mJy bins for the flux histograms, we re- 
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Figure 5. Left: The histogram of flux density values in the AzTEC/GOODS-S map. The dashed curve shows the distribution of 
flux values averaged over 100 noise realisations for this field, and is Gaussian distributed about zero. The solid curve shoves the flux 
distribution averaged over fully simulated maps, populated according to the best-fit Schechter function model to the GOODS-S data. 
This demonstrates how a fluctuation analysis is used to determine the best model to the GOODS-S data. Right: The grey-scale shows 
the likelihood values of S" and N^^jy from the fluctuation analysis, and the cross indicates the best-fit parameters to the data. The 
inner and outer contours indicate the 68.3% and 95.5% confidence regions, respectively. The error bars represent marginalized 68.3% 
confidence intervals on each parameter. 



strict the data-model comparison to bins with ^ 10 pix- 
els on average (flux densities between —2.8 and 5.5 mJy). 
The resulting best-fit parameters are {S'jN^rnjy) = 
(1.30l!|l:25mJy, leOljgmJy'Meg"^). The flux distribution 
for this best-fit model is shown as the solid curve in the left 
panel of Figure [5] The likelihood values for the S' — Nsmjy 
parameter space are shown in the right panel of Figure O 
with the best-fit parameters indicated by the cross. Due to 
the strong bin-to-bin correlations, it is not possible to de- 
termine the errors on the best-fit parameters analytically. 
Instead, we determine the errors statistically though sim- 
ulation by generating 600 fully simulated maps populated 
assuming the best-fit Schechter function model to the real 
GOODS-S map (including Poisson deviations), and then 
performing the same fluctuation analysis on these simulated 
maps. The distribution of best-flt parameters from these 
simulated maps are used to determine the 68.3% and 95.5% 
confidence intervals (contours in Figure [5]). The errors given 
for S' and A'smjy above (and shown as error bars in Fig- 
ure[Sl) represent the marginalized 68.3% confldence intervals 
on each parameter. 

The ability to recover the input number counts deter- 
mined from the fully simulated maps verifles the reliability 
of this fluctuation analysis method. The best-fit model to 
the actual GOODS-S data is listed in the first row of Ta- 
ble [2l and the differential number counts from this fluctu- 
ation analysis is shown in Figure [6] as the thick solid curve 
(best-flt) and dark shaded region (68.3% confldence inter- 
val). These results are also shown in Figures [7| and [8] for a 
comparison with other surveys. 

While potentially providing tight constraints on the 
source counts of the SMG population, this fluctuation anal- 



Table 2. The best-fit parameters for models to the 
AzTEC/GOODS-S number counts. The method used is listed 
in the first column: "P{d)" = fluctuation analysis, and "Bayes" 
= Bayesian method. The errors on the best-fit parameters rep- 
resent marginalised 68.3% confidence intervals. When an error is 
not listed, the parameter was fixed to the given value. 



Method 


Model 


S' 
(mJy) 


(mjy-ideg"2) 


a 


P{d) 
Bayes 
P(d) 
Bayes 


Eqn[3] 
Eqn[3] 
Eqn|4] 
Eqn|4] 


i.JU_o 25 
1 47+0.40 
J^-*'-0.25 


160l^^ 
107+?° 


-2.0 
-2.0 

"^-0.11 
o OC-+0.25 
■^■'^''-0.15 



ysis relies strongly on the accuracy of our noise realisations 
and the point source kernel, and is sensitive to the assumed 
number counts model. To demonstrate this, we carry out 
a fluctuation analysis on the GOODS-S data assuming in- 
stead a single power-law model for the number counts, which 
is given by 



dN 
d5 



iV: 



3m J y 



S 



3mJy 



a + l 



(4) 



We chose this functional form to make both parameters eas- 
ily comparable to those from the Schechter function model, 
since N^mjy gives the differential number counts at 3mJy, 
and a represents the power-law dependence at low flux den- 
sities. The best-flt parameters are listed in the third row 
of Table (2] and the best-flt model and 68.3% confldence 
interval are shown as the dashed curve and light shaded re- 
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Figure 6. The differential number counts for the 
AzTEC/GOODS-S field using various methods and models. 
The solid curve and dark shaded region indicate the best-fit 
model and 68.3% confidence interval from the fluctuation anal- 
ysis, assuming a Schechter function model. The dashed curve 
and light shaded region indicate the best-fit model and 68.3% 
confidence interval from a fluctuation analysis assuming a single 
power-law model. The squares show the Bayesian-extracted 
number counts and their 68.3% confidence intervals determined 
in Section 16.21 The slanted (horizontal) hatching shows the 
68.3% confidence interval for a fit to the Bayesian-extracted 
counts, assuming a Schechter function (power-law) model for the 
population. The horizontal dashed line shows the survey limit, 
where the number counts will Poisson deviate to OmJy~^deg~^ 
31.7% of the time given the area of the GOODS-S survey. 

gion in Figures [6l [71 and [S] The best-fit values of A^smjy 
for the two difTerent models are only marginally consistent, 
and the best-fit power-law index (q + 1 = —2.7) is signifi- 
cantly steeper than that assumed for the Schechter function 
model (of -I- 1 = —1.0). Comparing the two models in Fig- 
ure [51 they disagree significantly at S ~ 2mJy, where sup- 
posedly the tightest constraints can be placed on the number 
counts. One must thus exercise caution when interpreting re- 
sults from this fluctuation analysis, as it inherently assumes 
that the model provides a good representation of the source 
counts. 

6.2 Bayesian Estimation 

We can minimise the dependence of the analysis on our 
underlying model of the number counts by utilising a 
modified form of th e semi -Bayesian method introduced by 
ICoppin et al. I l|2005l . [ 20061 ). This semi-Bayesian method is 
now widely used to extract the number counts from sub- 
mm/mm surveys because of its ability to handle various 
survey biases, and it has been extensiv ely tested and val- 
idated using previous AzTEC data-sets (jPerera et al.|[20o3 : 
lAustermann et all [20091 , [20101 ). Since this method is de- 
scribed in detail in the aforementioned papers we only briefly 
summarise the steps here. 



The raw source counts from sub-mm/mm surveys suf- 
fer from three main biases: the "flux boosting" effect de- 
scribed in Section 14.31 contamination from false positives, 
and incompleteness. To account for the first two effects, we 
generate posterior flux distributions (PFDs), p{Si\Sm, o-m) 
(where 5"; is the intrinsic flux density of the source, Sm is the 
measured flux density, and am is the error on the measured 
flux density) , for each source candidate assuming some prior 
model for the SMG number counts. These PFDs are then 
randomly sampled (with replacement) to determine intrin- 
sic flux densities for the sources, and these fluxes are binned 
to calculate differential and cumulative number counts. This 
process is repeated 20, 000 times to adequately sample the 
number counts distribution. We also include sample vari- 
ance by Poisson deviating the number of sources sampled in 
each of the 20, 000 iterations. Since the PFD for each source 
candidate includes a non-negligible probability that the in- 
trinsic flux density is Si < OmJy, this procedure inherently 
accounts for false positives. 

To extract source counts from the AzTEC/GOODS-S 
map, we use the best-fit Schechter function model deter- 
mined from the fiuctuation analysis as the prior distribu- 
tion to generate PFDs for the source candidates. We sample 
all source candidates where the probability of the source 
having a flux density less than zero is p{Si < 0|5'm,cr„i) ^ 
0.05, which corresponds to S/N > 2.8 (actually depends 
on both Sm and am) fo r a total of 54 source candidates. 
[Austermann et al] l|2010l ) tested various limiting thresholds 
for p{Si < 0\Sm,am) using several AzTEC data-sets and 
found that any variations in the resulting number counts 
are much smaller than the formal 68.3% uncertainties. They 
found that for accurate PFDs, the p{Si < 0\Sm,am) Hmit- 
ing threshold does not greatly affect the resulting number 
counts (provided source confusion is not an issue) , and that 
using a higher limiting threshold supplies more information 
(due to increased survey completeness at faint flux densities) 
without introducing signiflcant biases. We have verifled that 
the PFDs for GOODS-S source candidates with S/N ^ 2.8 
are accurate to better than 1% at Si ^ 0.5 mJy f ollow- 
ing t he simulations described in [Austermann et al] ([2009I . 
[2010l ). further justifying the use of a p{Si < 0\Sm,am) ^ 
0.05 limiting threshold in this analysis. We note that while 
[Austermann et al] l|2010l ) use a limiting threshold of p{Si < 
[ Sm , a m ) ^ 0.2, we choose a more conservative limit to 
avoid potential systemat ics arising from confusio n effects, 
which were negligible for [Austermann et alj l[2010l ). 

The raw number counts must also be corrected for 
incompleteness. In previous implement ations of the semi- 
Bayesian method on AzTEC data-sets ([Perera et al.[ [2008[ : 
[Austermann et al.[|2009l , [2010l ), survey completeness was es- 
timated by the recovery rate of synthetic point sources with 
known intrinsic flux densities inserted into pure noise reali- 
sations one at a time. However, we have shown in Section[5]2] 
that this method overestimates the completeness for the 
AzTEC/GOODS-S held, since (unlike the AzTEC/ JCMT 
surveys) confusion noise is significant for our data. Using 
a pure noise completeness estimate would consequently un- 
derestimate the number counts in this field. We instead esti- 
mate the survey completeness through fully simulated maps 
as described in Section 15.21 where the simulated maps are 
populated according to the number counts distribution given 
by the assumed prior. To be consistent with our catalogue 
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selection from the real GOODS-S map, an input source in 
these simulations is considered to be recovered only if its cor- 
responding output source passes the limiting threshold test 
of p{Si < 0\S,n,am) ^ 0.05. Since we are using an "ideal" 
prior determined from the fluctuation analysis of this field, 
we are confident that this provides a good completeness es- 
timate for the correction of the raw number counts. 

The 1.1 mm differential number counts for the GOODS- 
S field determined from the Bayesian method are shown in 
Figures [5] and [7| (filled squares), and the cumulative num- 
ber counts are shown in Figure [S] The counts are also 
listed in Table O The number counts are calculated from 
the mean number of sources in each flux bin (with bin-size 
= 1 mjy) over the 20, 000 iterations, and the errors represent 
the 68.3% confidence intervals calculated from the distribu- 
tion in the counts across those iterations. For the differential 
number counts, the flux densities in Table [3] are the effec- 
tive bin centers weighted by the assumed prior. The number 
counts from this method are highly correlated since they 
are estimated by averaging over many realizations of the 
number counts bootstrapped off the same source catalogue. 
The covariance matricies for the differential and cumulative 
number counts are listed is Tables |4] and O respectively. 

The Bayesian-extracted number counts agree within the 
la errors with the best-fit Schechter function model from the 
fluctuation analysis. However, the number counts in the two 
lowest flux density bins (0.5 — 1.5 mJy and 1.5 — 2.5 mJy) 
are low compared to the best-fit curve from the fluctua- 
tion analysis, while the number counts in the highest bin 
(5.5 — 6.5 mJy) are high. This may arise from a systematic 
bias in the Bayesian-extraction method due to source blend- 
ing, since this technique does not account for the possibility 
that an individually detected source is the summed flux den- 
sity of two or more galaxies. This would indeed result in the 
apparent bias seen here, as the number counts would be 
overestimated at high flux densities and underestimated at 
low flux densities. We have checked for this potential bias in 
the Bayesian method by running this analysis on source cat- 
alogues extracted from 600 fully simulated maps populated 
according to the best-fit model from the fiuctuation analy- 
sis. Since in this case we know the exact form of the source 
counts input into each map, we can search for such effects in 
the output number counts. We find that the output differ- 
ential number counts for these simulated maps do indicate 
that this bias due to source blending is present in the data: 
for the 5.5 — 6.5 mJy flux density bin, the extracted number 

Table 3. The differential and cumulative number counts for the 
AzTEC/GOODS-S field, calculated using the Bayesian method 
described in Section 16.21 The errors indicate 68.3% confidence 
intervals. 



Flux Density dN/dS Flux Density N{> S) 

(mJy) (mJy-Meg-2) (mJy) (deg-^) 
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0.50 
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3.92 
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4.92 
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Table 4. The covariance matrix for the differential number counts 
for the AzTEC/GOODS-S field. The units are in mjy-^deg-*. 
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counts are higher than the input number counts for 65% of 
the simulated maps, while for the 0.5 — 1.5 mJy bin, the ex- 
tracted counts are lower than the input counts for 60% of 
the simulated maps. However, at all flux densities, the ex- 
tracted number counts agree with the input number counts 
within their la {2a) errors at least 86% (96%) of the time, 
so this bias is small compared to the formal Poisson errors. 

To verify that the Bayesian-extracted number counts 
are insensitive to the assumed prior, we reran this pro- 
cedure on the AzTEC/GOODS-S map using a prior dis- 
tribution that is consistent with a fit to the number 
counts from the AzTEC/SHADES survey; (5", iVsmjy, a) = 
(1.11 mJy, 153 mJy"Meg"2, -2.0). We find that the ex- 
tracted numbers counts using these two different priors agree 
within 4% at flux densities ^1.5 mJy. For the lowest flux bin 
at 0.5 — 1.5 mJy, the results agree within 19%, i.e. well within 
the formal la error. This demonstrates that the results from 
this technique are robust given a reasonable assumption for 
the prior number counts distribution. For this reason, we 
can fit these number counts to various models. 

For a given model, we fit to each of the 20, 000 bootstrap 
iterations separately as the fl ux density bins for a given it- 
eration are uncorrelated (see lAustermann et al.l 2010), and 
we use the distribution of best-flt parameters to determine 
the most likely values and their confidence intervals. The 
results of a fit to the GOODS-S differential number counts 
assuming a Schechter function model (Equation O and a 
power-law model (Equation [Jjl are given in rows 2 and 4 of 
Table [2] and are shown in Figure |6] as the hatched regions, 
which indicate the 68.3% confidence intervals. The fits to 
the Bayesian-extracted counts are in good agreement with 
the results from the fluctuation analysis for a given model. 
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To compare the two models, we compute the metric given 
by 

= {d — m) w (d — m)''" (5) 

where d is the row-vector containing the differential number 
counts from the Bayesian method, w is the corresponding 
weight matrix calculated from the inverse of the covariance 
matrix, and m is the model number counts evaluated at the 
same flux density bins as d. For the best-fit Schechter func- 
tion to the Bayesian-extracted number counts, ~ 0.84. 
This model provides a better fit to the data than the single 
power-law model, for which ~ 5.4. Note that since the 
errors are not Gaussian-distributed, this metric is not ex- 
pected to be follow the x^-distribution; we use this metric 
simply to compare the relative goodness-of-fit for these two 
models. 



6.3 Comparison with Results from SCUBA 
Lensing Cluster Surveys 

The AzTEC/GOODS-S survey currently provides the only 
constraints on the 1.1mm number counts at Si.imm < 
1 mjy. For comparison, none of the existing blank-field 
SCUBA surveys constrain the 850 ^.m counts to compar- 
ative ly faint fiux densiti e s (Ss50iiui < 2mJ y): the deepest 
(e.g. iHughes et all Il998l : lEales et alTl2000l ') are too small 
in area to provide statistically significan t samples of faint 
SMGs, while the large-area surveys (e.g. iBorvs et al.ll2003l : 
ICoppin et al.l i2006) do not reach sufficient depths to probe 
such faint sources. The most significant constraints on the 
Sgso/jm ^ 2mJy number counts come from SCUBA lensing 
cluster surveys, where massive foreground clusters are used 
to probe faint background SMGs via gravitational lensing 
(e.g. Small ct al, 2002; Cowicjyt _aL, .2002 ; Chapm an et al.l 
120021 : iKnudsen et all \20(M . \200^ . This is a powerful tech- 
nique for detecting intrinsically faint high-redshift SMGs 
and for estimating their source counts at low flux densi- 
ties, as the foreground clusters magnify both the flux of the 
background sources (by factors of typically 2 — 3) and the 
area in the source plane, effectively decreasing the survey 
confusion-limit that hinders the sensitivity of blank-field ob- 
servations. With new constraints on the 5*1. imm = 0.5 mJy 
number counts from the AzTEC/GOODS-S survey, the faint 
end of the number counts determined from lensing cluster 
surveys can be compared to the results from a blank-field. 
Given the small sample size from lensing cluster surveys, 
this comparison is limited to the cumulative number counts. 

The number counts from four separate SCUBA lensing 
cluster surveys are shown in Figure |8] To compare these 
with our 1.1mm results, we must account for the differ- 
ence in observing wavelengths. We assume here that SCUBA 
and AzTEC are sampling the same underlying population 
of sub-mm/mm-bright galaxies, and that the difference in 
the observed number counts can be described by a simple 
flux scaling: R — S'ssoMm/'S'i.imm. In principle, surveys at 
1 mm may preferentially select sources at higher redshifts - 
or sources with colder dust temperatures - than surveys at 
850 [im. While there is some evidence that 1 mm surveys se- 
lect on average higher redshift galaxies than 850 \im surveys 
llEales et all i200a : lYounger et all |2007| . l2009l : iGreve et all 
|2008| ). other studies suggest that there is no signiflcant 



difference between the two population s (jCreve et al.ll20o3 : 
llvison et al1l2005l : iBertoldi et al.ll2007l ). 

A recent source-to-source compari son of SCUBA and 
AzTEC sources in the GOODS-N held (jChapin et al.ll2009l ) 
reveals that while the redshift distribution of 1.1 mm sources 
in that field peaks at a higher redshift than that of 
850 [J.m sources {z — 2.7 versus z = 2.0), the popula- 
tion is consistent with an average flux scaling with R = 
1.8. Alternatively, the 850 \im numbe r counts determined 
from the SCUBA/SHADES survey (ICoppin et al.1 I2006D 
and the 1.1mm number counts from the GOODS-N field 
(Pcrcra ct al. 2008) are consistent assuming a flux scaling of 
7? = 2.1 ±0.2. These estimates are equivalent to the expected 
flux ratio of a 2 = 2.5 galaxy whose SED can be modelled as 
a single component modifled black-body with T^j = 30 K and 
emissivity index /3 = 1.5. As this model is consistent with the 
expected SEDs of local galaxies observed with the Infrared 
Astronomical Satellite (IR AS) and SCUBA (jPunne et all 
I2OOOI : iDunne fc Ealesll200ll) a s well as the measured SEDs 
of several SMGs llChapman et al. 20051: iKovacs et al.|[2006l : 
iPope et al.ll2006l : ICoppin et al.ll2008 l. we start by adopting 
a simple scaling factor of 7? = 1.8 for the purposes of com- 
paring the number counts from the SCUBA lensing cluster 
surveys to the AzTEC/GOODS-S number counts. 

The scaled number counts from the SCUBA lensing 
cluster surveys agree only marginally with the number 
counts from the AzTEC/GOODS-S survey and are system- 
atically higher at all flux densities (see Figure [S]). In fact, 
the lensing cluster survey number counts are systematically 
higher than those from all AzTEC blank-field surveys. We 
first examine the possibility that the flux sc aling factor of 
R = 1 .8 derived ft-om the GOODS-N field bv lChapin et all 
(|2009l) is not universal by deriving the value of R that min- 
imises the residuals between the SCUBA lensing cluster 
number counts and the best-fit model to the 1.1 mm number 
counts from the AzTEC/GOODS-S field. The best-fit values 
and 68.3% confidence errors for R determined from each of 
the four lensing cluster data-sets shown in Figure|8]are given 
in Table [S] We find R ^ 3.3 for all lensing cluster surveys 
considered here, regardless of which model (Schechter func- 
tion or power-law) we consider. Such high values of R are 
inconsistent with observed values for individual bright SMGs 
(when both 850 \im and 1.1 mm measurements are available) 
as well as predicted values assuming SEDs and redshifts typ- 
ical of this population. We note however that a high value 
of R — 2.5 ± 0.1 has been estimated from a similar scaling 
of the 850 [im and 1.1 mm num ber counts from the SCU BA 
and AzTEC/SHADES surveys l|Austermann et al.|[201oh . If 
1.1 mm surveys are tracing a population of galaxies with ei- 
ther higher redshifts or colder dust temperatures than those 
detected at 850 |i.m, we would expect this flux ratio to be 
lower, since these two bands probe closer to the dust peak. 
The high value of R measured from scaling the SHADES 
number counts is likely due to systematics caused by differ- 
ences in the number counts analyses or calibration, and we 
consider this an upper limit to the true flux ratio. Estimates 
of R from scaling the 850 [im lensing cluster number counts 
to the best- fit model to the AzTEC/GOODS-S data-set rep- 
resent > 2cr deviations from this upper limit oi R = 2.5 from 
the SHADES results, implying that even greater systematics 
exist in the lensing cluster number counts. Evidence for sig- 
nificant systematics in the derivation of the lensing cluster 
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Figure 7. The differential number counts from the AzTE C/GOODS-S surve y {black squ ares), compared with th ose determined 
from AzTEC surveys of other fields, including: GOODS-N llPerera et al.ll2008^ . COSMOS l|Austermann et al.llioogh . and SHADES 
dAustermann et al.H201ol) . The error bars represent 68.3% confidence intervals on Bayesian-extracted counts. The solid (dashed) curve 
and dark (light) shaded region indicate the best-fit Schechter function (power-law) model and 68.3% confidence region from a fluctuation 
analysis of GOODS-S (Section [6TTJ. The horizontal dashed line shows the survey limit, where the number counts will Poisson deviate to 
OmJy-ldeg-2 31.7% of the time given t he area of the GOOD S-S survey. The 95% confidence interval from a fluctuation analysis of the 
1.1mm Bolocam Lockman Hole survey l lMalonev et al.ll2005h assuming a single power-law model is shown by the hatched region. The 
filled circles indicate the differential number counts determined from a Bayesian analysis of the 870 |J.m LABOCA survey of the ECDF-S, 
scaled to 1.1mm assuming a flux ratio of SgTOfim/Si imm = 2.0. The dotted curve shows the best-fit Schechter function model to the 
LABOCA/ECDF-S data from a P{d) analysis scaled to 1.1mm (confidence region not available). See the on-line journal for a colour 
version of this Figure. 



number counts is also given by the fact that the cumulative 
number counts at S'sso/jm > 2 mJy from lensing cluster sur- 
veys are also higher t han those from large-area blank-field 
SCUBA surveys (e.g. iBorvs et al.ll2003l : IWebb et al.l l2003l : 



IScott et al.]|2006l : ICoppin et al.ll200g ) 

There are several issues in the number counts extraction 
from lensing cluster surveys which may systematically bias 
the number counts high. Errors in the detailed mass models 
of the lensing clusters are an obvious candidate for the ob- 
served bias; however, any errors in estimating magnification 
factors for source flux densities are co mpensated by equiv- 
alent errors in the source plane area (|Knudsen et al.ll2006l . 
I2OO8I ). Contamination by foreground cluster members is an- 
other concern, but these have been identified and excluded 
in these lensing cluster analyses. All of the lensing cluster 
surveys considered here include sources detected with low 
significance {S/N ^ 3.0), and thus could include a signifi- 
cant number of false positives which would bias the number 



counts high. However, we believe that the dominate sys- 
tematic that causes the discrepancies between our number 
counts and those derived from these lensing cluster surveys 
is due to the fact that non e of these analyses, excepting 
that of iKnudsen et af] (|2008l ). ac count for the fi u x boo sting 
effect described in Section 14.31 iKnudsen et al.l l|2008l ) ap- 
ply a correction for fiux boosting only for the sources de- 
tected in their de epest maps (^'ss o um 5. 2mJy). If we fit for 
R using only the IKnudsen et al.l l|2008l ) number counts at 
<S'850Vim ^ 2mJy (row 5 in Table |6]), the best-fit value agrees 
with our conservative upper limit of i? = 2.5 within the 
1(7 errors. However, this effect arises from the steep num- 
ber counts distribution for SMGs: it is more likely that a 
source detected in a fiux-limited survey is an intrinsically 
faint source (numerous) boosted high by noise than an in- 
trinsically bright source (rare) . Flux boostin g is thus promi- 
nent in all maps, and we believe that the IKnudsen et al.l 
(2008i' ) number counts have not been fully corrected for this 
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Figure 8. The cumulative number counts from the AzTEC/GOODS-S survey and AzTEC surveys of other fields. See the caption of 
Figure [7] for a full description of the symbols. The 1.2 mm number counts from [Creve et al are shown as empty upside-down 

triangles with 95% confidence error bars. For comparison, the number counts from SCUBA 850 [xm tensing cluster surveys are shown 
with smaller symbol s. All 850 |xm count s have been scaled t o 1.1 mm assuming a si mple flux ratio scaling o f S'ssOMm/S'i.imm = 1.8. The 
number counts from I Cowie et aT || |2002|'| .[ Small et al. I ll2002l' ), IChapman et al.1 j2002l'l . and lGreve et al" have not been corrected for 

flux boosting are therefore are likely biased to the right. See the on-line journal for a colour version of this Figure. 



Table 6. The 850 H-m to 1.1mm flux ratio estimated by scal- 
ing the number counts from SCUBA tensing cluster surveys to 
the best-fit model to the AzTEC/GOODS-S number counts. The 
columns are: 1) reference for the tensing cluster survey (last row 
is the results for all four surveys combined); 2) best-fit flux ra- 
tio R = SssOMm/Sl.lmm and 68.3% confidence itervat assuming 
the Schechter function model for the GOODS-S number counts 
(Equation Is} ; and 3) best-fit fiux ratio and 68.3% confidence in- 
terval assuming the power-law model for the GOODS-S number 
counts (Equation [4]| . 



Reference R (Schechter) R (Power-law) 



Cowie et al. (2002) 
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effect. Still, we note that if SMGs cluster strongly on small 
angular scales, the apparent bias in lensing cluster surveys 



may simply be the result of cosmic variance. The total num- 
bers of lensing clusters analysed is still small, and the errors 
on the cumulative counts for the lensing cluster surveys are 
dominated by Poisson noise due to the limited survey areas 
(s; 40arcmin^ in the source plane) and small sample sizes. 



6.4 Comparison with Other 1.1mm Surveys 

The differential and cumulative number counts from other 
1.1 mm blank- field surveys are shown in Figures [7] and [8] for 
com parison. Thes e include the AzTEC s urveys of G O ODS- 
N fPerera et al.' '2008'), COSMOS (IScott et al.' l2008l: 
Aust crmann et al. 2009), and SHADES ( Austcrman n et al.l 
201(1 ) taken on the JCMT. The SHADES data-set consists 
of two separate regions of sky - the Lockman Hole-East and 
the Subaru/XMM-Newton Deep Field - and covers a total 
area of 0.67 deg^, making it the largest blank-field survey at 
1.1 mm published to date. The 95% confidence interval from 
a fiuctuation analysis of the 1.1 mm Bolocam Lockman Hole 
surve y (la = 1.4 mJv: iLaurent et al. I l2005l : iMalonev et al.l 
l2005h . where a single power-law model is assumed, is also 
shown. In Figure[8l we show the 1.2 mm cumulative number 
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counts derived from the M AMBQ surveys of the Lockman 
Hole and ELAIS N2 fields (|Greve et al.ll2004l ). 

The best-fit single power-law model for the 1.1mm 
number counts determined from our fiuctuation analy- 
sis of the AzTEC/GOODS-S map is in good agreement 
with that determined from the fluctuation analysis of the 
Bolocam Lockman Hole survey. However, given that the 
AzTEC/GOODS-S survey is 2 - 3 times deeper than the 
Bolocam LH survey with roughly the same area, we are 
sce ptical of the ext r emely tight confidence intervals reported 
by iMalonev et al. I (|2005l ). The single power-law model has 
been shown to poorly describe the SM G number counts 
over a wide range in flux density (e.g. | Scott et al.l l2006l : 
ICoppin et al.ll2006l : I Austermann et al.ll201ol v with the best- 
fit slope likely a compromise between the steeper slope at 
high flux densities and the shallower slope at low fluxes. The 
resulting model inflexibility is a likely contributor to the 
unde rrepresented e rrors. While the 1.2 mm number counts 
ft-om [g reve et al] (|2004l ) are broadly consistent with the 
number counts determined from several AzTEC surveys, 
we note that this analysis does not fully account for the 
substantial flux boosting effect on the measured source flux 
densities. 

The general agreement among the 1.1mm blank-field 
surveys taken with AzTEC is quite good. The identical, 
well-tested methods used to reduce and analyse these data- 
sets minimises any systematic differences among the re- 
sults, so this agreement is naturally expected. For the same 
reason, any subtle differences seen among these surveys is 
significant. The number counts from the GOODS-S field 
are consistent with those from the SHADES survey, which 
currently provides the tightest constraints on the blank- 
field 1.1mm number counts over this flux density range. 
In contrast, the GOODS-N and COSMOS flelds appear 
somewhat over-dense by as much as ~ 50%; we see a sig- 
nificant amount of variation in the number counts on ~ 
100 arcmin^ scales. Since galaxies are organized within struc- 
tures with significaiit pow er on comoving scales < 100 Mpc 
(e.g. ISpringel et ai]|2005l ). variations in the observed num- 
ber counts is expected for surveys covering < Ideg^Q A 
quantitative study of the degree of field-to-field variations 
observed among all blank-field surveys taken with AzTEC 
on the JCMT and the ASTE will be presented in a future 
paper (Scott et al. in prep.). 



7 COMPARISON WITH THE LABOCA 
SURVEY OF THE ECDF-S AT 870 

A 0.35 deg^ region that covers the Extended Chan- 
dra Deep Field South (ECDF-S) has been imaged at 
870 n m with the LABO CA camera on the APEX tele- 
scope (|Weifi et al.ll2009l ). The LABOCA ECDF-S Submil- 
limetre Survey (LESS) therefore encompasses the entire re- 
gion mapped by AzTEC at 1.1mm and reaches a uniform 
complementary depth of la ~ 1.2mJybeam~^ at 870 |J.m. 
A detailed comparison involving a joint analysis of the two 
data-sets by a combined group from the two surveys is in 



^ 10' = 5.1 Mpc at 2: = 2 in the concordance model of cosmology. 



progress. In this Section we carry out a preliminary compar- 
ison of the sub-mm/mm sources and their number counts in 
order to check for general consistency between these two 
data-sets. 



7.1 AzTEC- and LABOCA-Detected Sources 

In Figure [9l we show the AzTEC 1.1mm map with the lo- 
cations of the AzTEC and LABOCA sources. Of the 126 
870 iJ.m source candidates detected with S/N ^ 3.7 in the 
LABOCA catalogue (where five false detections are ex- 
pected), 20 are located within the 50% uniform coverage 
region of the AzTEC survey. Using a conservative search 
radius of 30" (roughly the FWHM resolution of both sur- 
veys), we find that 16/20 LABOCA sources are also de- 
tected by AzTEC with S/N ^ 3.5. Turning this around, 
16/41 AzTEC source candidates are detected by LABOCA 
with S/N ^ 3.7. Additionally, aU seven of the S/N ^ 4.5 
AzTEC source candidates located in the extended (i.e. nois- 
ier) region of our map are also detected by LABOCA. 
AzTEC so u rces with probable 870 iJ.m counterparts from 
IWeifi et all l|2009l l are noted in Table □ AzTEC/GS42 and 
AzTEC/GS43 each coincide with pairs of LABOCA source 
candidates separated by < 25". 

Verifying that the degree of overlap between the two 
source catalogues is consistent with expectations given the 
completeness of the surveys is not straightforward given 




Figure 9. The AzTEC/GOODS-S flux density map with the 
positions of AzTEC and LABOCA source candidates indicated. 
To distinguish between higher- and lower-significance detections, 
AzTEC sources detected with S/N > 5.0 are shown by larger 
circles, while AzTEC sources detected with S/N < 5.0 are shown 
by smaller circles. LABOCA sources detected with S/N ^ 5.0 
are indicated by crosses, and LABOCA sources detected with 
S/N < 5.0 are indicated by pluses. The contour outlines the 50% 
uniform coverage region for the AzTEC/GOODS-S survey. 
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the different methods used to extract sources and quan- 
tify the map properties for these different surveys. This 
wiU be fully addressed in the joint analysis of both data- 
sets. The primary reasons for the lower fraction of AzTEC 
sources detected by LABOCA are 1) we use a lower S/N 
threshold than that used for the LABOCA catalogue, 
and 2) the AzTEC map (scaled in ffux density assuming 
•S'sTO/jm/'S'i.imm ~ 2) is slightly deeper than the LABOCA 
map. We stress that we do not expect a large fraction of the 
AzTEC-identified sources to be random noise peaks. 

For 21 sources that are listed in both the LABOCA 
and AzTEC source catalogues, we compute the 870 |i.m to 
1.1mm flux ratio from the deboosted fluxes (we exclude 
AzTEC/GS42 and AzTEC/GS43 since each are resolved 
into two LABOCA sources). The mean flux ratio of this sam- 
ple is (S'sTOMin/S 1 . 1mm ) = 2.0 ± 0.6. Scaling from 870 ^.m to 
850 [im assuming a spectral index of /3 = 1.5 (5*^ oc ly^'^), 
this corresponds to a flux ratio of (S'gsp^jm/iS'i. imm) = 2.1 ± 
0.6, c onsistent with previous estimates from IChapin et al.l 
l|2009l ) and from scaling the 850 [im and 1.1mm num- 
ber counts derived from other surveys (|Perera et al.|[200^ : 
lAustermann et al] l2010t ) . The observed flux ratios of the 
AzTEC/LABOCA sources are also consistent with expecta- 
tions for cold {Td = 25 — 50K), dusty galaxies at 2 ^ 1 — 3. 
From this direct source comparison there appears to be no 
signiflcant systematic differences (e.g. calibration) between 
these two data-sets. Note however that since this estimate of 
the average 870 [im to 1.1 mm flux ratio is calculated solely 
from LABOCA-detected AzTEC sources, it may be biased 
high if there is a significant population of very high-redshift 
(or very cold) dusty galaxies that are faint at 870 [im. 

7.2 870 vim and 1.1 mm Number Counts 

IWeifi et al.l (|2009l ) carry out a P{d) analysis on this data-set 
assuming four different parametrised models. Since all of 
the best-fit models give nearly the same results for the mag- 
nitude and shape of the number counts, we consider only 
the Schechter function parametri sation for the purposes of 
comparmg the IWeifi et al. I l|2009l ) results to ours. We note 
that since the LESS field covers a much larger area than 
the AzTEC/GOODS-S field, we might expect some differ- 
ences in the source counts from these two data-sets owing 
to cosmic variance. The best-fit Schechter function to the 
LABOCA/ECDF-S data, scaled to 1.1mm assuming a fiux 
ratio of (Sgsojim/S'i.imm) ~ 2.0 (as measured from the direct 
source comparison in Section [7. If) , is shown in Figures [7] and 
[1 The number counts from the LABOCA/ ECDF-S field 
derived using the semi-Bayesian method of ICoppin et al.l 
l|2006l) are also included in these figureslf] 

Considering first the results from the fiuctuation analy- 
ses, we see fairly good agreement between the number counts 
derived from the LESS and AzTEC/GOODS-S surveys until 
5*1. imm ^ 1.0 mjy, at which point the LESS number counts 
rise more steeply. This is largely du e to the diff erences in 
the analyses: while we fix a = -2.0, IWeifi et all |2009) fit 
for this parameters and find a best-fit value of a = —3.7 
(as defined by Equation As explained in Section 16.11 

^ The Bayesian-derived number counts from LESS shown in Fig- 
ures [7] and |8] were provided by A. Weifi. 



we choose to fix a to a value consistent with that mea- 
sured by previous AzTEC surveys because the GOODS-S 
data cannot well-constrain a 3-parameter model. It is likely 
that the larger LABOCA/EC DF-S survey is able to put use- 
ful constraints on a; however IWeifi et al. do not in- 
clude error estimates on their best-fit parameters, so it is 
not possible to quantitatively access the degree to which 
our results are consistent. We note that the value of a de- 
termined from the LESS fiuctuation analysis is well-matched 
to our best-fit spectral index from a fiuctuation analysis of 
AzTEC/GOODS-S assuming a single-power law model. 

Comparing the Bayesian-extracted number counts from 
the AzTEC/GOODS-S and LESS surveys, we find that they 
are in excellent agreement for 5*1. imm ^ 3mJy; however 
at Si.imm ^ 3mJy, the Bayesian-extracted number counts 
from the LABOCA data are significantly lower than our esti- 
mates. ,Weifi_etj|L, (. 2009 ) noted a similar difference between 
their P{d) and Bayesian-estimated counts and speculated 
that this discrepancy at low fluxes arises from source blend- 
ing, since two or more faint sources blended together by the 
large beam would be counted as a single bright source and 
result in an underestimate of the n umber count s of SM Gs at 
low flux densities. For this reason IWeifi et all (|2009l ) place 
more confldence in their P{d) results. However we have 
shown through simulations that the P{d) method can be 
highly model dependent. As described in Section 16.21 our 
simulations indicate that we largely account for confusion 
effects in our Bayesian-extracted number counts and that 
they are not significantly biased low (high) at faint (bright) 
fiux densities due to source blending in the catalogue. We 
thus place more confide nce in our Bayesian-derived results, 
bmce IWeifi et al.l l|2009l ) do not present similar tests of ei- 
ther of their number counts extraction algorithms, it is not 
possible to know the underlying cause of the discrepancies 
between their two different methods. 

IWeifi et all (|2009l ') report that the 870 ^m number 
counts from LESS are a factor of ~ 2 lower than those 
from the 850 ^m SCUBA/SHADES survey at flux densi- 
ties SgTOfim > 3 mJy, and note that this seems consis- 
tent with the reported under-densities in the ECDF-S of 
other high-redshift galaxy populations. On the contrary, we 
find that the 1.1 mm number counts from AzTEC/GOODS- 
S are completely consistent with those measured in 
the AzTEC/SHADES survey. Since sub-mm/mm-selected 
galaxies potentially span a much larger volume of space 
than high-redshift galaxies selected by other photometric 
criterion, we would not necessarily expect to see an under- 
density of SMGs in this field. As the number counts derived 
from the LESS, AzTEC/GOODS-S, and AzTEC/SHADES 
fields are largely consistent assuming a reasonable scaling 
for the differences in observed wavelength, it seems un- 
likely that cosmic variance alone can account for the fac- 
tor of ~ 2 difference in the number counts derived from 
LESS and SCUBA/SHADES. Considering also the physi- 
cally unrealistic high flux-scaling factor required to bring the 
SCUBA/SHADES and AzTEC/SHADES number counts 
into agreement {R = 2.5), one possibility is that this dis- 
crepancy arises from systematics in the SCUBA/SHADES 
ca libration and /or number counts extraction as discussed 
bv lAustermann et al.l (|2010l ): this will be properly addressed 
via a source to source comparison of AzTEC and SCUBA 
sources in the SHADES fields (Negrello et al., in prep.). 
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8 CONTRIBUTION OF DIFFERENT GALAXY 
POPULATIONS TO THE COSMIC 
INFRARED BACKGROUND AT 1.1 MM 

The FIRAS instrument on the Cosmic Background Ex- 
plorer (CORE) detected a uniform cosmic in frared back- 
ground (CIRB), which p eaks around 200 (i.m (jPueet et alJ 
Il996l : iFixsen et allll998h . The total power in the infrared 
background is 1 — 2.6 times larger than that of the opti- 
cal/UV background, suggesting that the majority of radia- 
tion produced by young stars and AGNs over t he cosmic his- 
tory is absorbed and reprocessed by dust (e.g. ICispert et all 
I2OOOI ). If dusty starburst galaxies in the early universe ac- 
count for much of the CIRB, the AzTEC/GOODS-S survey 
is expected to resolve a significant fraction of the CIRB. 

Summing the deboosted 1.1mm flux densities of the 
41 S/N ^ 3.5 source candidates in the AzTEC/GOODS-S 
map, we measure an integrated flux of 1.5±0.2 Jy deg"'^ over 
the 0.075 deg^ field. Comparing this to the total energy den- 
sity in the CIRB a,t 1.1 m m of 18 - 24 Jy deg"^ ()Puget et all 
ll996l : lFixsen et al.lll998h . we have resolved only 6-8% of the 
CIRB into individual galaxies. However, if we instead inte- 
grate the best-fit Schechter function model to the GOODS-S 
number counts down to mjy, we estimate that we have sta- 
tistically detected 6.3 Jydeg"^, or 26 — 35%, of the CIRB at 
1.1mm with our AzTEC/GOODS-S survey. The fact that 
we do not resolve 100% of the CIRB through this integration 
implies that simply extrapolating the Schechter function 
model that best fits our data at 0.5 mJy < 5i.imm < 6.5 mJy 
to much lower flux densities significantly underestimates the 
number counts at the faint end. On the other hand, inte- 
grating the best-fit single power-law model to the GOODS- 
S number counts down to SLimm ~ 0.04 mJy would ac- 
count for 100% of the CIRB at 1.1mm; however, we know 
this model to be a poor description of the number counts 
of bright SMGs as well as physically unfeasible since it ap- 
proaches infinity as 5'i.imm goes to zero. 

We next use a stacking analysis to estimate the frac- 
tion of the CIRB at 1.1mm that is resolved by the en- 
tire 1.4 GHz radio population. Stacking at the locations of 
N = 222 radio sources in this field (up slightly from 219 
radio sources used in the stacking analysis in Section 14.21 
since we have shifted our AzTEC map to correct the as- 
trometry), we calculate an average 1.1mm fiux density of 
5'i.imm, radio ~ 660 ± 78 ^J.3y ■ Assumiug that each of the ra- 
dio sources distributed over an area of A = 0.075 deg^ has 
a fiux density of 5i,imm, radio, the radio population has a to- 
tal integrated fiux of N ■ 5'i.imm,radioM = 1-9 Jydeg"^, or 
8 - 11% of the CIRB at 1.1 mm ' 

We finally estimate the contribution to the CIRB 
at 1.1mm from MIPS 24 jxm-selected sources using a 
similar stacking analysis. The average fiux density from 
1185 24 jxm sources distributed over a 0.068 deg'^ area is 
Si. imm,24/^m — 290 =b 26 /^Jy. The total integrated fiux from 
24 Jim sources at 1.1mm is 5.0 Jydeg"^, or 21 — 28% of the 
total CIRB. This is similar to the fraction of the CIRB 
at 850 vim resolved by 2 4 [oa sources (29 — 37%) in the 
SCUBA/GOODS-N field (|Wang et al.ll2006l ). In contrast, a 
stacking analysis of 24 yim sources with the BLAST maps 
of GOODS-S at 250, 350, and 500 ^m suggests that the 
full intensity of the CIRB at these sho rter wavelengths is 
resolved by sources selected at 24 [im l|Devlin et al.1 l2009l : 



iMarsden et"ai] |2009| ). This possibly demonstrates the ex- 
istence of a significant population of higher-redshift (z > 
3) dust-obscured galaxies that are (statistically speaking) 
missed by current A < 500 \mi surveys, but account for 
~ 2/3 of the CIRB at longer wavelengths. 



9 CONCLUSIONS 

We imaged a 270 arcmin^ field towards the GOODS-S re- 
gion to a confusion-limited depth of la « 0.6 mJy using the 
AzTEC camera on the ASTE, making this the deepest sur- 
vey carried out to date at 1.1 mm. We detect 41 SMG candi- 
dates with S/N ^ 3.5, where roughly two are expected to be 
false positives arising from noise peaks. This survey is 50% 
complete at 2.1 mJy and 95% complete at 3.5 mJy. We have 
demonstrated that the presence of confusion noise has signif- 
icant consequences for the properties of the map (and hence 
the number counts) and must be considered when accessing 
the survey completeness and expected number of false de- 
tections in the source catalogue. From realistic simulations 
of the SMG population in this field, we estimate that f» 26% 
of the source candidates identified in the AzTEC/GOODS- 
S map are actually two or more mm-bright galaxies with 
comparable fiux densities blended together due to the low 
angular resolution of the ASTE beam. 

We have used two very different methods to estimate 
the SMG number counts in this field: a fiuctuation analy- 
sis where we model the distribution of fiux density values 
in the map, and a semi-Bayesian technique where the num- 
ber counts are determined from sampling the PFDs from 
the catalogue of SMGs. We have demonstrated that both 
methods are able to retrieve the correct number counts dis- 
tribution from fully simulated data-sets. Furthermore, the 
best-fit number counts to the GOODS-S field using these 
different methods are consistent. We note that while our 
fluctuation analysis can provide good constraints on the 
number counts, the results depend strongly on the assumed 
model for the number counts distribution. The depth and 
large survey area of the AzTEC/GOODS-S map have re- 
sulted in the tightest constraints to date on the SMG num- 
ber counts at 5'i.imm = 0.5 mJy. Comparing our cumu- 
lative number counts to those from several SCUBA lens- 
ing cluster surveys at 850 (xm, the lensing cluster number 
counts appear to be biased high assuming reasonable values 
for the flux-scaling from 850 (xm to 1.1mm. These results 
are consistent with those from large-area blank-field sur- 
veys at 850 \im with SCUBA, where the number counts at 
Sssofii-n 2mJy from lensing cluster surveys are systemati- 
cally higher; this is most likely due to fiux boosting effects 
not being fully treated in the number counts extraction from 
the lensing cluster surveys. We find that the number counts 
from the AzTEC/GOODS-S field are consistent with those 
from the 0.67 deg^ AzTEC/SHADES survey. 

Comparing our source catalogue to the 870 vim source 
list from the LABOCA/ECDF-S survey, 16/20 (80%) 
LABOCA sources located within the 50% uniform cover- 
age region of the AzTEC/GOODS-S survey are also de- 
tected with S/N ^ 3.5 in the AzTEC map; however since 
the AzTEC map is slightly deeper, only 16/41 (39%) of the 
AzTEC sources in this field are listed in the LABOCA source 
catalogue. In contrast to the results from the LABOCA sur- 
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vey, we find no apparent underdensity of SMGs in tfiis fieid 
compared to previous submm/mm surveys. For tlie SMGs 
tliat are fisted in both catafogues, tfie mean 870 \im to 
1.1mm flux ratio is 2.0 ± 0.6 and is consistent with our 
expectations given typical dust SEDs and redshifts for these 
gafaxies. Scaling the 870 vim number counts to 1.1mm, we 
find fairly good agreement between the counts for Si.imm ^ 

1.0 mjy (> 3.0 mjy) extracted using a fluctuation analysis 
(Bayesian method) on each of the two surveys. We think 
that the discrepancies between the AzTEC and LABOCA 
source counts at fainter flux densities are due to differences 
in the assumed models used for the fluctuation analyses, and 
additional biases in t he Bayesian metho d for the LABOCA 
analysis as noted by IWeifi et all (|2009l '). A combined anal- 
ysis of these two confusion-limited surveys at 870 |J.m and 

1.1 mm is in progress. 

We resolve only 6 — 8% of the CIRB at 1.1 mm into in- 
dividual mm-bright galaxies. While the 24 |i.m population 
can account for the full energy density in the CIRB at 
250 — 500 (im, we estimate that 24 fim sources statistically 
detect only 21 — 28% of the CIRB at 1.1 mm, demonstrating 
that a population of faint dust-obscured galaxies at « > 3 
that are largely missed at shorter wavelengths potentially 
contribute significantly to the total energy density in the 
CIRB at 1.1mm. 
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